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Abstract 

The existence of large Vn — v T mixing suggests the likelihood of large smuon-stau mixing in su- 
persymmetric models, leading to /i and t number violation. In addition to interesting signatures 
in slepton and neutralino production and decay, this will lead to rare r decays, such as r — * jj/y. 
Recently, it has been pointed out that the r — > 3/i branching ratio could be substantial in the 
large tan/? region of parameter space, due to an induced fi — t— Higgs vertex. In this paper, 
another signature, t — > [ir\ is considered. In the large tan (3 region, it is shown that the branching 
ratio of r — * fir] is 8.4 times the branching ratio of r — > 3/i, independent of any unknown param- 
eters, and it will thus give the most stringent bound on Higgs-mediated lepton flavor violation, 
and may provide its first signature. In the other regions of parameter space, where r — > fij is 
the most prominent decay, the branching ratio for r — > \ir\ is always substantially lower. 

The flavor physics of quarks and leptons is one of the most prominent mysteries in particle physics. 
The most surprising development in flavor physics in the past decade has been the observation |Q] of 
very large mixing between muon and tau neutrinos. The mixing will, at some level, lead to mixing 
in the charged lepton sector, giving violations of muon and tau lepton number conservation. 

In the Standard Model supplemented with right-handed neutrinos, such violation is generally 
small 0| . However, a much bigger effect will occur in supersymmetric models, by inducing mixing 
between the scalar muon and the scalar tau. In the most general supersymmetric standard model, 
even if neutrino mixing were not present, one would have large mixing between all scalar leptons and 
scalar quarks. This would lead to very large flavor (quark and lepton) changing neutral currents, 
which are not observed. It is thus generally assumed that the squark and slepton masses are equal 
at the unification scale; an assumption that occurs naturally in many models, such as supergravity 
or gauge-mediated supersymmetry breaking. 

Even if the masses of the smuon and stau are equal at the unification scale, without mixing, 
the presence of non-diagonal neutrino mass terms (either from different Dirac mass terms or non- 
diagonal right-handed neutrino mass terms) will affect the masses through renormalization group 
running, and will generate mixing termsQ. In this paper, we will consider only mixing between 
left-handed smuons and staus; in most models, such mixing is the largest. It should be kept in mind 
that solar neutrino oscillation experiments HJ also indicate large mixing between muon and electron 
neutrinos, and in models with inverted hierarchies, there could be substantial mixing between left- 
handed smuons and selectrons, although the very strong bounds on fx — > e^f will constrain these 
effects. 



The mixing is characterized by the parameter £23 = M| 3 /m 2 , where M23 is the off-diagonal term 
in the slepton mass matrix and m is the slepton mass scale. The value of S23 is extremely model- 
dependent, of course, but in models in which the mixing arises entirely through renormalization 
group running H |6|, 0, its value is typically between 0.1 and 1. 

The most studied tau-number and muon-number violating process is r — > ^[O]. Many of these 
works consider various models for <!>23- Normalizing the rate to the current boundj8| |9| 

In the constrained MSSM, where the parameter space is restricted to a manageable dimensionality, 
this process dominates in the low mjft regionQ. 

In addition, one can consider tau-number and muon-number violation in production and decay 
of sleptons, neutralinos and charginos[||, [ji], |H| . As an example, the process \2 — > Xi^ T where xi,2 
are neutralinos, can be searched for at the LHC. In the constrained MSSM, this process dominates^ 
in the 777,1/2 > ma region of parameter-space, and is thus complementary to r — > fi'y. 

Recently, Babu and KoldaQ pointed out that r — > 3fi was a promising signature in models 
with a large value of tan/3. Earlier Jl3[, they had noticed that squark mixing would induce a flavor 
non-diagonal quark-quark-Higgs Yukawa coupling, and had examined the consequences for rare B- 
decays. The same mechanism, however, will also induce a /j — r — Higgs vertex, and thus directly 
to r —* 3/j, through tree-level Higgs exchange (either the h, H, or A). The branching ratio, for a 
reasonable choice of mass parameters, is 

, . ~ /tan/3\ 6 /100GeV\ 4 

where m,A is the pseudoscalar mass. This result is very insensitive to the SUSY spectrum, with the 
exception that it can increase by up to a factor of four for large /j. This branching ratio will be 
accessible at B-factories in the near future. 

In this paper, another signature of lepton-number violation is discussed: r — ► firj. It will be 
shown that, in the large tan (3 region discussed in the previous paragraph, the branching ratio is 
much higher than r — > 3/U, and may be a much more sensitive test of Higgs-induced lepton flavor 
violation. 

In the Babu-Kolda model(5|, the lepton flavor violating interaction is given by 

C LFV = (2G 2 F )^^l(r Rf , L )[H° + iA°] + h.c. (2) 

where H° and A are the heavier scalar and the pseudoscalar, respectively, and we have chosen 
the generally preferred region of parameter-space in which sin(a — /3) ~ 1. Here, K32 is given in 
Ref. H and depends on loop integrals (but is relatively insensitive to SUSY parameters). For 
SUSY parameters \i = M\ = 771,2 = m \ = fn^ui Mr = 10 14 GeV and the off-diagonal Dirac 
neutrino coupling equal to the top quark Yukawa coupling (as expected in 50(10) models), one has 
K32 = 4 x 10" 4 . 

With this interaction, one can have a r convert into a \i and a virtual H° or yl , which then 
converts into a f^ + ^" pair. This gives the branching ratio mentioned above. However, one could 
equally well have an A convert into a strange quark pair, which then becomes an 77, as shown in 




Figure 1: Diagram leading to r 
and Kolda. 



[if]. The small circle is the lepton-flavor violating vertex of Babu 



Figure 1, giving r — * fir]. This will have both advantages and disadvantages. The two body phase 
space is a major advantage, and the extra color factor and slightly bigger coupling (since m s > m„) 
are also advantages. The disadvantage is in converting the strange quarks into an r\. 

A general discussion of r — * fir] can be found in Ref. [14]. The relevant matrix element is 



< 0|s75s|?7 >= — V6F 1 



v m u + md + 4m s 



(3) 



Using this matrix element, it is straightforward to calculate the branching ratio. If one divides by 
the r —* 3/i branching ratio, the unknown parameters all cancel, and the result is (neglecting the 
muon mass) 



T(r -» 3/i) 



m,, 



mi 



(4) 



With F° ~ 



150 MeV[15], this ratio is 8.4, giving a branching ratio of 



BR(r -» firj) = (0.84 x 10" 



tan (5 
60 



100 GeVV 



m A 



One can get this result approximately without doing a calculation. Imagine that final state 
interactions are turned off and that the strange quarks propagate as free particles. One expects the 
ratio of t — > fiss to r — ► 3/i to have a factor of 3 for color and a factor of {m s /m^) 2 for the Yukawa 
coupling. The cross diagram in the muon case turns out to lower the rate by a factor of 3/2, so the 
overall ratio is ~ 10. Since the only two body decays would be un and uri', and the latter is 
suppressed much more by phase space, the fir) rate should dominate this process. 

Since the experimental bound (T^| on the branching ratio for r — > \ir\ is 9.6 x 10~ 6 , and thatQ 
for t — ► 3/i is 1.9 x 10~ 6 , it is clear that r — ► fir] puts stronger constraints on the model. In order 
to reach the interesting region of parameter space (tan/3 ~ 60 and uia ~ 100 GeV), the bound on 
t — > 3/i would need to be improved by a factor of 20, whereas the bound on r — > fir] would need to 
be improved by a factor of 10. 

Could these improvements be made? Both the r — > 3fi and r — * fir] bounds are based on the 
CLEO-II sample of 4.7 fb _1 , and the CLEO experiment has now accumulated a total of 5 times 
that luminosity (which would give a total of 24 million tau pairs). So in the absence of backgrounds, 




Figure 2: Box diagram leading to r — > /it/. A similar diagram will lead to t — > /i7r. The are 
neutralinos. 



the current bounds could improve by a factor of five. The efficiency for r — » 3/i is listed as 15%; 
the efficiency for r — > /i?? is about 3%, when one includes the fact that they only search for the 
77 channel for the r\ (thus the factor of five difference in the current bounds). Including the 
three-pion decay, or increasing the fiducial area for finding photons, could improve that efficiency 
substantially. Over the next years, BABAR and BELLE will reach 500 fb _1 , which could easily 
reach the needed sensitivity, depending on the point at which they become background-limited. 
Even if tan (3 is somewhat smaller, or wia larger, the necessary sensitivity could possibly be reached 
at LHC, SuperKEKB or a tau-charm factory. Note that the r — > 3fi decay could still be dominant 
in the small region of parameter space in which Mjj « Ma- 

Are there any other processes that can give r — > fir]? One can have the box diagram of Figure 
2, which will also yield r decays into fi plus other mesons, including the tt, p and 4>. If we take the 
special case in which the neutralinos are pure photino, the rate for r — > fin is given by 



r( r - fin) _ 32™F 2 m§(/ 2 + if) 



81m2m?Mf (x) 



(5) 



where x = m^/mj, is the photino mass, mr is the average slepton mass, M%{x) is given in Ref. 
||, and the integrals are 



16vr 2 



dx 



(x, \x 2 ) 



(x — l) 2 (x — a)(x — b)(x — c) 



(6) 



Here, a, b, c are the ratios of the squark,smuon and stau masses to the photino mass. In giving this 
expression, we have used the fact that m^/(m u + md) is numerically close to m T . In the case of 
r — > fir], there will be a suppression of a factor of 4 from the s-quark charge, an increase of a factor 
of {F~t I F w ) 2 ~ 3.0 from the decay contants, and the coefficient of the I 2 term will be decreased by 
a factor of (3m^/8m 2 ) relative to m 2 , or a factor of about 3. This ratio has been evaluated for the 
entire SUSY parameter space, assuming sparticle masses in the range of 5 — 1000 GeV[^8[, and is 
always less than 10~ 2 . As a result, r — ► /i+ meson arising from the box diagram is negligible. 

The existence of large — v T mixing implies mixing between the left-handed smuon and stau. 
While one could look for this mixing directly in neutralino/slepton interactions, one can also look 
at r decays. The decay r — * fvy is one signature, however Babu and Kolda have noted the mixing 
will also lead, especially in the large tan/3 region, to r — > 3]i. In this paper, it has been pointed out 
that t — > fit] will also occur in this Higgs-mediated model, with a branching ratio 8.4 times bigger, 



and is thus more sensitive. In other models, where r — > ^7 is the main signature, the r — ► fir] rate 
is substantially smaller. 
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